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I. INTRODUCTION 



One of the most important and intriguing recent developments in theoretical physics is the AdS/CFT correspondence. 
An example of holographic gauge/string theory duality [li] asserts the complete physical equivalence between particular 
string vacua and ordinary quantum field theories. In certain regimes it relates strongly coupled field theories to weakly 
curved string theory, which in the leading order reduces to supergravity. A very interesting class of applications involves 
systems at finite temperature. The holographic duality relates thermodynamic notions in quantum field theory to 
black hole mechanics in the bulk description. Even more fascinating are non-equilibrium phenomena which cannot 
be described by thermodynamics or static black hole solutions. There is very strong motivation for such studies on 
both sides of AdS/CFT duality. 
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While there are several frameworks that provide fundamental explanations for the entropy of a stationary black hole, 
the entropy of dynamical black holes is an even more difficult problem. This issue, in view of AdS/CFT duality, is 
directly connected with the physics of non-equilibrium processes in strongly coupled quantum field theory. Progress 
in this area will have important bearing on problems in black hole as well as particle physics. 

The past few years have seen tremendous activity applying gauge/gravity duality to study gauge theories beyond 
the perturbation series. One obvious application area where such a tool is sorely needed is the investigation of non- 
perturbative dynamical QCD plasma. The case of A/" = 4 super Yang-Mills (SYM) theory has been the focus of 
attention in this context, since its holographic representation in terms of string theory on anti-de Sitter spacetimes 
is the best understood example of gauge/gravity duality^. Furthermore, it appears that the significant differences 
between this theory and real-world QCD do not play a major role in a particular range of temperatures relevant to the 
heavy ion experiments currently in progress at RHIC and soon to start at the LHC. Thus the studies of Yang-Mills 
plasma at finite temperature using tools provided by string theory (such as are currently of 

great practical importance and attract a lot of attention from the heavy ion community. 

It is now fairly well established that there is a regime where the dynamics of quark-gluon plasma is well described by 
relativistic hydrodynamics. The hydrodynamic description itself is well understood despite the fact that the underlying 
quantum field theory is strongly coupled. It is based on the idea of an expansion in the number of gradients, much in 
the spirit of effective field theory. One spectacular application of AdS/CFT in recent years was to show that in the 
case of A/" = 4 SYM a hydrodynamic description can be derived by considering time-dependent black-brane solutions 



in AdS spacetime (particularly useful review is [lO|). Further applications to quark-gluon plasma require a better 



understanding of the duality at finite temperature and in a dynamical setting. 

A byproduct of the hydrodynamic construction is the notion of entropy current. Recent advances in relativistic 
conformal hydrodynamics show that a phenomenologically introduced hydrodynamic entropy current consturucted 



order by order in the gradient expansion involves an ambiguity [8|, At second order a 4-parameter family of 



currents^ was identified in [s^. Subsequently some important considerations on this topic have appeared in the 
literature fill, which suggest that it may be possible to reduce this freedom to just one parameter. It seems however 
that there is a real ambiguity in the hydrodynamic entropy current and in the light of AdS / CFT duality it is natural 
to ask if one can understand its origins on the gravity side. One of the goals of the present study is to address this 
issue. 

At the root of this question lies the identification of the area increase theorems in general relativity with the second 
law of thermodynamics. Thus if one hopes to understand this ambiguity on the gravity side, the first step is to 
carefully examine and understand how horizon areas increase. In doing this one may build on the experience gained 
in the study of standard black holes in asymptotically flat four-dimensional spacetime. Generalizing these results to 
five-dimensional anti-de Sitter spacetime is not difficult, but the deep questions encountered earlier remain. 

The most "conservative" definition of entropy identifies it with the area of a spatial section of the event horizon. 
This notion has the drawback of being a teleological (and thus global) concept. In the context of AdS/CFT duality 



^ Many of the results hold also for sectors of other large Nc strongly coupled gauge theories which have gravity dual described by the 
Einstein-Hilbert action with negative cosmological constant. This is another hint pointing towards possible (limited) applicability of 
AdS/CFT results to the real world. 

^ In general there might be a 5-parameter ambiguity, but one of these parameters is connected with parity-odd effects. This article 
considers only the parity-even case. 
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this leads to acausality in the field theory 1^, |l3[ . In classical general relativity, the problems raised by the global 
character of the event horizon have lead to alternative, quasilocal notions of black holes. These include trapping |l4| . 
isolated [3], and dynamical horizons. These programmes are closely related to each other and motivated by 
historical ideas about trapped surfaces and apparent horizons fisl . 

These quasilocal horizons have a problem of their own, namely non-uniqueness: unlike event horizons, in dynamical 
black-hole spacetimes, there are many possible time-evolved apparent (or trapping or dynamical) horizons. In the 
case of apparent horizons, each foliation of the spacetime will give rise to a different time-evolved horizon [l^ while 
trapping/dynamical horizons are also subject to deformations (see, for example, [3] 20|). Though this is generally 



thought to be a bad thing, in the context of AdS/CFT it is natural to suspect that the ambiguity of the hydrodynamic 
entropy current may be related to ambiguities of this type. 

In this paper it is shown that the non-uniqueness of apparent horizons is not the source of the ambiguity in the 
entropy current. Instead, the already existing uncertainty as to whether the event or apparent horizon determines 
the entropy is embraced and a more general approach is advocated which mimics the phenomenological construction 
of the boundary entropy current. In the spirit of the membrane paradigm j2l| we consider "horizons" that are made 
up from families of (not necessarily trapped) codimension-two surfaces that satisfy properties such as area increase, 
asymptoting to the correct equilibrium limits, and being "almost" apparent horizons. In general such families are hard 
to deal with, but the situation becomes much more manageable close to equilibrium. We consider spacetimes in which 
a small parameter can be used to formulate a perturbative expansion around equilibrium (this does not necessarily 
imply a hydrodynamic description 1^). Such a formalism - slowly evolving geometry - has been developed in the 
context of regular four-dimensional black holes 
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123,1231 and it is adapted to the case of black brane mechanics 
as well as generalized to the case of not-quite-apparent horizons. This formalism also makes it possible to see the 
conditions under which the first law of thermodynamics can be formulated in a dynamical setting. 

It is conjectured here that the ambiguity in the above definition of black brane entropy corresponds to the known 
ambiguity of the hydrodynamic entropy current in the appropriate regime on the gravity side. Verifying this claim 
in the general case is beyond the scope of this paper, which explores this question in a special, highly symmetric 
case - the gravity dual to Bjorken (boost-invariant) fiow. Besides being tractable, this geometry has the important 
feature that there is a unique apparent horizon consistent with symmetries of boundary flow. This means that the 
ambiguity in the hydrodynamic entropy current in case of Bjorken flow cannot be interpreted as a consequence of 
slicing dependence. 

Boost-invariant flow is one of the simplest, yet phenomenologically interesting examples of boundary dynamics. It is 
a one-dimensional expansion mimicking the dynamics of plasma created in the heavy ion collision with an additional 
assumption of boost-invariance along the collision axis. This is based on Bjorken's observation 2J] that multiplicity 



spectra when expressed in proper time and rapidity variables are approximately independent of rapidity in the mid- 
rapidity region. While this approximation may be somewhat rough, it provides a dramatic simplification. After 
taking into account all the symmetries in this problem, it turns out that all the physical quantities depend on single 
variable - the proper time. In the pioneering work Janik and Peschanski considered the gravity dual to boost- 
invariant fiow in the regime of large proper time and showed that regularity of the bulk geometry forces perfect fluid 
hydrodynamics on the boundary. Subsequent developments included calculating subleading corrections to the large 
jroper time geometry, which correspond to dissipative terms in the boundary energy-momentum tensor Q, 3| (see 



25l | for a useful review). Understanding the gravity dual to Bjorken flow within the framework of fiuid-gravity duality 
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27| makes it possible to address issues of black brane mechanics in this case^ 



One can also look at this paper from another perspective. In understanding any mathematical formalism, it is very 
useful to have concrete examples with which to work. There are very few exact (or even perturbative) dynamical 
black holes solutions on which to test ideas such as those about slowly evolving horizons. Most of the already known 
examples were examined in {2^. These included spherically symmetric spacetimes where the expansion was driven 
by matter flows and perturbations of Schwarzschild where the expansion was driven by shears (gravitational waves) 
- though unfortunately in this case the order of the perturbation was such that one could see the shears driving the 
expansion but not the expansion itself. Then, from this point of view, the current work presents a new non-numerical 
example of a near-equilibrium horizon and it is the first one for which a shear driven expansion can actually be 
directly observed. This allows equations to be checked in previously inaccessible regimes and also helps to build a 
better understanding of the formalism. 

The structure of the paper is the following. Section two reviews the various definitions of black hole and then develops 



black hole mechanics in AdS„+i spacetime, generalizing [20|, |22| and keeping things as general as possible with a 



view to further applications of these results. The third section gives a brief introduction to both field theory and 
gravity aspects of the boost-invariant flow and the fourth section applies methods developed in section two to the 
boost-invariant situation. The fifth section contains the analysis of the freedom in the definition of entropy on both 
sides of the AdS/CFT correspondence in the case of Bjorken fiow and shows that there is precise match. The final 
section offers conclusions and signals some directions for further research. 



II. BLACK BRANE MECHANICS 



A. A brief review of black holes and horizons 



We begin with a review of three ways of thinking about black holes. Since our main interest in this paper are the 
dynamically evolving black branes dual to Bjorken flows, we will focus in particular on dynamical black holes. While 
the various definitions of black holes generally agree for stationary black holes, they diverge away from equilibrium. 
For example, they identify different surfaces in spacetime as the boundary of the black hole region and these surfaces 
have different surface areas. This will be significant in later sections when we compare the thermodynamics on the 
CFT side with the corresponding black hole mechanics on the gravity side. Then, as noted in the introduction, 
different surface areas would suggest different values of entropy. 



Event horizons: causally defined black holes 



It is well known that causally defined black holes and event horizons are not locally defined. Paraphrasing, a classically- 
defined black hole [3| is a region of spacetime from which nothing can ever escape. At first glance this seems like 
a perfectly reasonable definition, however a little consideration quickly turns up problems and these arise from the 



It is critical to have a smooth description of the horizon region. 
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FIG. 1: A schematic demonstrating the non-local nature of event horizon evolution for a spherically symmetric spacetime with 
the angular dimensions suppressed. Horizontal location measures the radius of the associated spherical shell while time is 
(roughly) vertical. The shaded gray region represents infalling null dust. 



false event horizon 




true event horizon 



concepts of "ever" and "escape" . To identify a black hole region one must essentially sit at infinity and wait forever 
to make sure that all escaping signals are identified and further that those that initially look like they might escape 
really do make it to infinity. Equivalently (but more rigorously) the black hole region is the complement of the causal 
past of (future null infinity). The boundary of the black hole region is the congruence of null geodesies known as 
the event horizon. 

To better understand this definition and its associated peculiarities consider Figure [T] which represents the evolution 
of an initially isolated Schwarzschild black hole which is later irradiated by an infalling shell of null dust (see 23] for 
more discussion of the Vaidya spacetimes used to generate this example). In the diagram the inward null direction 
is horizontal while the outward null directions are tangent to the light gray dashed lines. Thus in general the future 
of any point is "up-and-to-the-left" with the light cones pointing more and more towards the singularity as one 
approaches it. 

In this case one simply tracks the evolution of radial null geodesies to find the event horizon ~ if they fall into the 
singularity they are inside the black hole but if they are still heading outwards after the shell of matter passes then 
they are deemed to escape to infinity. Intially, without any knowledge of the future arrival of the shell of matter, one 
would guess that the heavily dashed line at r = 2Afo (Afo the mass parameter in the initial Schwarzschild spacetime) 
would be the event horizon: null geodesies inside it (labelled A in the diagram) clearly head inwards towards the 
singularity while r — 2Mo itself keeps a constant area. However, with the arrival of the matter shell one finds that 
this false event horizon also terminates at r = as do apparently escaping geodesies like B. It then turns out that the 
true event horizon was a congruence of geodesies that initially appeared to be escaping but was later "almost caught" 
by the increased gravitational field arriving with the shell. After the passage of the shell these become stationary 
r = 2(Mo + AM), where Mo + AM is the new Schwarzschild mass. 

This example nicely demonstrates the consequences of a causal structure definition of black holes. While the evolution 
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of any congruence of null geodesies is certainly causal, the identification of the set corresponding to the event horizon 
depends on events in the far future. As a result non-omniscient observers cannot precisely locate them. Further even 
if identified they will evolve in non-intuitive ways: in the example, the arrival of matter didn't cause the event horizon 
to expand but instead curtailed an expansion that began earlier (in a sense in anticipation of the mass increase) . 

Apparent horizons: geometrically defined black holes 

Given these peculiarities, it is often more useful to turn to alternative definitions of black holes and their boundaries 
which leave aside the causal structure of spacetime and instead focus on the strong gravitational fields characterized 
by the existence of trapped surfaces. For regular four-dimensional astrophysical black holes, trapped surfaces are 
closed and spacelike two-surfaces which have the property that all families of null geodesies that intersect them 
orthogonally must converge into the future. To understand this intuitively, consider a transparent spherical shell that 
is covered with light bulbs and sitting in empty space. Then if the bulbs are quickly turned on and then off again, 
two spherical light fronts will be generated - an outwards moving one that expands in area and an inwards moving 
one that contracts. By contrast, if the shell is transported so that it lies inside a Schwarzschild black hole, concentric 
with the horizon and enclosing the singularity, then both light fronts will fall towards the centre of the black hole and 
contract in area; again consider Figure [T] where both outward and inward falling congruences contract in area. This 
is the canonical example of a trapped surface. 

More mathematically, if and n° are respectively the outward and inward pointing null normals to a two-surface S 
then one can write 

6'(£) < and 6'(„) < , (1) 

where 0(£) and are the expansions of the null normals and are geometrically analogous to (traces of) extrinsic 
curvatures 

Ow = r'yJb and 0(„) = r'Van, (2) 

with 

qab = 9ab + Lnb + riaib (3) 

being the induced (spacelike) metric on the two-dimensional surface. Alternatively given outward and inward con- 
gruences of null geodesies which have tangents i"" and n"- on S one can show that 

V^6'(£) = ^Av^and v/^6'(„) = i£„v/^, (4) 

where y/q is the area element on S and C indicates is the Lie derivative operator. Then it is clear that the sign of the 
expansion determines whether the congruence is expanding or contracting. 

More generally, given the energy conditions the mere existence of a trapped surface in an asymptotically fiat spacetime 
implies both: 1) a singularity somewhere in its centre (future) [17] and 2) that it is necessarily contained in a causal 
black hole and so event horizon [18]. Indeed, for the standard Kerr family of (stationary) black hole solutions, the 
set of all points contained on some trapped surface coincides exactly with the black hole region. Thus it is not 
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FIG. 2: An "instant" Ef along with some of its trapped surfaces (small black circles), the associated trapped region (dark gray) 
and the apparent horizon (thick dashed line). 




unreasonable to consider the existence of trapped surfaces as being the key characterizing feature of a black hole 
region and this is the basis of the alternative definitions of black holes. 

The original such definition was the apparent horizon. This begins with the foliation of a spacetime into spacelike 
hypersurfaces - essentially instants in time. Then at a given instant, the trapped region is the union of all the trapped 



surfaces contained in the hypersurface and the boundary of that region is the apparent horizon. It can be shown 18| , 
that on the apparent horizon = and 9^^) < 0. More generally, any such surface satisfying these conditions is 
referred to as marginally trapped'^. 

This observation that apparent horizons are marginally trapped then motivates the various modern notions of quasilo- 
cal horizons such as trapping [l^ . isolated [3], and dynamical horizons [igI (or see review articles such as 29 1 and 



Though there are significant technicalities, the idea is that the identification of marginally trapped surfaces, 
which under arbitrarily small deformations become fully trapped, is sufficient to signal the presence of a black hole 
boundary - even without going through the process of foliating the spacetime and finding apparent horizons on each 



slice. In fact, this idea is so pervasive that in numerical relativity (31|, the term "apparent horizon" has been co-opted 
to refer to the outermost surface 6'(£) = on a give slice of spacetime. The study of these ideas is an active and 
developing area of research with a fairly complicated system of nomenclature but for the purposes of this paper, 
quasilocal black hole (or brane) horizons in a (n + l)-spacetime dimensions will be understood as n-dimensional 
hypersurfaces that are foliated by {n — l)-dimensional marginally trapped spacelike surfaces. The requirement that 
there be fully-trapped surfaces "just inside" the horizon can be mathematically written as: 

< , (5) 



where in this case n° is any extension of the null normal into a neighbourhood of the putative horizon [20| . In what 
follows we will usually adopt Hayward's nomenclature and refer to such structures as future outer trapping horizons 
(FOTH) [3]. Occasio nally however we will refer to time-evolved apparent horizons (which are examples of FOTHs) 
or dynamical horizons |16| which are almost equivalent. 

Mathematically the properties of these surfaces (such as existence, uniqueness, and evolution) may be studied using 
standard techniques from differential geometry. From a physical perspective however, the key advantages of these 
generalized apparent horizons include: 1) they are defined by the existence of strong gravitational fields, 2) as for 
fully trapped surfaces their existence is sufficient to imply the existence of singularities and event horizons, 3) they 



* Unfortunately there is some inconsistency in the literature about the usage of this phrase. That said, our usage here is probably the 
most common one — for more details on usage see [29| . 
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FIG. 3: A simulation similar to that of Figure [T] though this time two distinct shells fall into the black hole. Both the apparent 
and event horizons are plotted. 
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may be identified without reference to the far future (if one thinks of them as time-evolved apparent horizons), and 
4) their evolution is similarly local. The local nature of this evolution is demonstrated in Figure [H In that case, the 
expansion of the event horizon continues to occur in anticipation of the arrival of infalling matter, with the actual 
arrival of the matter slowing or ending that expansion. By contrast the apparent horizon evolves in the expected way 
in response to the infalling matter - it expands in and only in the presence of actual matter crossing the horizon. 

On the downside, it is well-known that quasilocal horizons are not uniquely defined. For classically defined apparent 
horizons this is easily be seen. Given a foliation of spacetime we can define a time-evolved apparent horizon A as the 
union of the apparent horizons on each surface. Then, it is clear that different foliations will sample different subsets 
of all the possible trapped surfaces. Thus, different foliations will define different A. In the extreme case it is known 
that certain slicings of Schwarzschild spacetime contain no trapped surfaces at all and so no apparent horizon 
We will return to this lack of uniqueness in later sections. 



32|. 



The membrane paradigm: a physical approach to black holes 



A third way of looking at black holes focuses not on causality or geometry but rather on how black holes interact with 
their environment. By definition event horizons cannot directly affect their surroundings (they are not in causal contact 
with any point outside themselves) and neither can apparent horizons (they are contained within event horizons). All 
that either can do is impose restrictions on the behaviour of surfaces "near" the horizons that are in causal contact 
with the outside. The membrane paradigm [21] studies the physics of n-dimensional timelike hypersurfaces "just 
outside" the event horizon. For astrophysical purposes, it turns out that one can view these as the time evolutions 
of (n — l)-dimensional spacelike viscous fluid surfaces which carry energy, angular momentum and entropy and can 
exchange these with their surroundings. The details of this formalism are not important here but it is important 
to understand that for some purposes surfaces besides the usual horizons can usefully be thought of as physically 
characterizing black holes. 
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FIG. 4: A schematic of an n-tube A with compact foliation surfaces S\ along with the outward and inward pointing null 
normals to those surfaces. V" is the future-pointing tangent to A that is simultaneously normal to the 5*^ . 




B. The geometry of n-tubes 

With the motivation of the last section in mind, the mathematics of n-tubes, a class of geometric structures that 
includes both FOTHs and event horizons as well as the time-like hypersurfaces of the membrane paradigm, will be 
examined here. In an (n -I- l)-dimensional spacetime n-tubes are n-dimensional surfaces which can be foliated by 
(n — l)-dimensional spacelike surfaces ^a- The term "tube" comes from (3 -|- l)-dimensions where the S\ for horizons 
are compact and generally diffeomorphic to S"^ (Figure |4]), however this nomenclature will be kept even for black 
branes where the S\ are diffeomorphic to M.^ (or in the case of AdSn+i with R"^^) and so certainly not compact. 

Foliated event horizons, time-evolved apparent horizons and membranes are clearly examples of such structures. Event 
horizons are n-tubes of null signature which have the correct causal properties as discussed in the previous section. 
Time-evolved apparent horizons are FOTHs whose S\ are the apparent horizons found in individual space-time slices. 
As will be clear in moment, these are either null (if isolated and in equilibrium) or spacelike (if dynamical and 
expanding). Finally in the membrane paradigm the A is a timelike surface and the S\ can be thought of as the 
evolving configurations of a fluid. 

To begin, consider the basic geometry of n-tubes and in particular focus on the spacelike S\. First, the co-dimension 
of the is two and the normal space has Minkowski signature so, as in the brief review, one can always find null 
normals £a and Ua which span that normal space. is taken to be outward-pointing (and so is tangent to A if it is 
null) while n° points inwards towards the singularity; again see Figure [H For convenience the normals are usually 
cross- normalized so that £ ■ n = —1 which leaves a single scaling degree of freedom in their definition 

r ^ /r and n" jn" (6) 

for any positive function /. 

The intrinsic geometry of the Sx is defined by the induced metric qat © while the extrinsic geometry is characterized 
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by the derivatives of the normal vectors in the directions tangent to S\. These inchide the extrinsic curvatures: 

fcifc = Q^fb^dd and fci;) = raq^VcUd , (7) 
which decompose into traces 9(^£) and 9(^„) as well as trace-free parts a'fjj and a'^^: 

Kb = (^^_^^ ^Wlab + <Jab andfc^^ " („_ ;^) %»)g°'' + ^afc ' (8) 

Physically these are respectively the expansions and shears of congruences of null curves which have tangents and 
n° as they intersect S\. 

The rest of the extrinsic geometry of the S\ is described by the connection on the normal bundle: 

[da^-qUcVbt. (9) 

The geometric information contained in is somewhat obscured by its gauge-dependence on the scaling /: 

£^ f£ and n ^ ju =^ iua ^ Cua + daln f , (10) 

where da is the (n — l)-dimensional gradient operator on S\. In the usual way for gauge potentials, the invariant 
information is contained in its curvature: in this case this is the curvature of the normal bundle 

Q.ab = daCOb - dbLUa ■ (H) 

These normal-bundle quantities will arise in the following discussions and physically are closely related to angular 
momentum (see for example the extensive discussion in [20|). That said, in the actual application of this work to our 
boost-invariant spacetimes, LUa will vanish thanks to the symmetry of the black-branes. 

The next natural step is to consider how the 5a fit together to form A. Let V° denote the evolution vector field 
tangent to A that maps leaves of the foliation into each other {CyX 1) and is normal to each of the S\. Then / 
can always be chosen so that 

V° = r - Cn" (12) 

for some function C. This ties the scaling of the null vectors to the foliation and the freedom is reduced to that in the 
foliation labelling: that is / = /(A). In particular, from ^TU\i it is straightforward to see that this fixes uja, removing 
its gauge dependence. 

Clearly thanks to the normalization £-n = — l,C>0<^Ais spacelike, while C = <^ A is null, and C < A is 
timelike. Further the "time" -rate of change of the area element^ is 

Thus, for apparent horizons with = and 0(„) < 0, C also characterizes the evolution of the horizons and is often 
referred to as the expansion parameter. Specifically in such cases 

C < <^ y/q is decreasing <^ V° is timelike, 

C = ^ ^/q is unchanging <^ V° is null, (14) 
C > 4=> \fq is increasing <^ V° is spacelike. 



^ Quotation marks are used around the word time since if V" is spacelike then this is a coordinate rather than physical notion of time. 
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In contrast, for event horizons C = but away from equilibrium > (thanks to the second law) and so the 
horizon can still expand. For the timelike surfaces of the membrane paradigm none of 0(„) or C vanish. Equation 
()13p will be very important in the following, also in cases when is not exactly zero. 



Another gauge quantity 

Kv = -V'^ni.Vat (15) 
describes how the null normal evolves up A. Under rescalings 

£ ^ f{\)i =^ ^ /«v + ^ . (16) 

In the case of an isolated horizon Ky — Ke, the familiar surface gravity that appears in the first law of black hole 
mechanics and it is not hard to show that a zeroth law of black hole mechanics holds: Kg is necessarily constant over 
the horizons [3] (or section Hi CI of this paper), though its exact value is fixed by the chosen normalization of the null 
vectors. 



Constraints on the geometry 



These geometric quantities are not all independent. Instead they are linked by a set of constraint equations that 
ensure that the S\ fit smoothly together to form an n-tube and further that that tube must be embeddable in a larger 
(n + l)-dimensional spacetime (on which the Einstein equations hold). We now consider some of these constraints. 

First, the "time" rate of change of the two-metric on the can be written in terms of the extrinsic curvature 

CvQab = qlqfCvqcd = 2(^11^ - Cfci^) ) = (0(,) - C6^n))qab + 2{af^ - Cai';' ) (17) 
and this is consistent with the already discussed (|13p if one takes the trace of both sides. 

More involved are the equations relating the derivatives of extrinsic curvature quantities to the rest of the geometry. 
Here we just list some of these results that will be used in our upcoming discussion of mechanics. Those interested in 
more details of the derivations of these results can consult Appendix El First, for 6'(£) 



and 



(fC + 2Cj''daC - C 



- + \\cJ(i)f + Gabtt 



(18) 



For e 



(n) 



n- 1 



in) 



\\uf + daC^'' - R/2 + Gabn'^t ~ 0(,)0(„) 



+c 



n — 1 



(19) 



(20) 
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and 

n — 1 

In these equations, Iplp = '^''li'a, ||cr*^^'"'' |P = ^ Ricci scalar on the (n — l)-surfaces, Gab = 

TZab ~ ^T^9ab is the Einstein tensor while Cabcd is the Weyl tensor. The only other equation that we will need is the 
"time" rate of change of oja- 



Gbcii' + Cn-) - CbcdeV^i^'n' 



CvCba = daKv - k[^^ - k[f [d^'C - [b''C] (22) 



(n-l) 



As an immediate application of equation Eq. (|18p. it can be used to (infinitesimally) quantify the lack of uniqueness 
of apparent horizons. Consider a particular marginally trapped surface So and fix a scaling of £"■. Then one can try 
to solve (1181) for C to find a direction in which So may be evolved while maintaining 0(f) =0. If a solution exists then 
it is unique, but different scalings of will give rise to different Cs and so different possible directions of evolution; 
in general. So will be a member of an infinite number of different marginally trapped n-tubes. That said, (|18p does 
strongly constrain the geometry of these possible evolutions from So- It can be shown that if one chooses a particular 
Ao, all other possible evolutions must lie partly in the causal past and partly in the causal future of A,,. Equivalently, 
all other A must "interweave" with the original one (the original proof of this result may be found in 19], but also 
see [2^ for discussion in a language closer to that used here). In particular if Aq is highly symmetric, then other 
A cannot share that symmetry. Thus, in highly symmetric spacetimes it is usually possible to select a preferred 
marginally trapped n-tube that shares those symmetries. 



C. Equilibrium states and the zeroth law 



We now begin an examination of how the geometric calculations of the last section give rise to the laws of black hole 
mechanics. First consider equilibrium states and the zeroth law. 

Equilibrium states are characterized by null n-tubes on which 0(f) = and —C^j^t^ is future-pointing and causal 
(a condition slightly weaker than the dominant energy condition). These are a type of isolated horizons known as 
non-expanding horizons [ist and as we now shall see include apparent horizons that are in (possibly temporary) 
equilibrium with their surroundings as well as event horizons in eternal equilibrium (such as the Kerr family of black 
holes) . 



First if A is null, then V" — 1°" and C — Q. Then (|T8|) reduces to the Raychaudhuri equation 

£f0(f) = tivO(t) - Ik^'^lP - Gab^" - el^/{n - 1) , (23) 



or with 0(f) ~ 0: 



lk*^^ir + Ga6«''-0. (24) 



On the geometry side it is clear that a\^^ = and so the intrinsic geometry of A is time invariant 

Ctqab = \k'^2 = if^W'iab + <J^ab = «■ (25) 
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Meanwhile, on the energy side, Raychaudhuri impHes that Gab^"'^^ = and so by the energy condition —G'^i^i'' = fil"' 
for some function /i: there are no matter flows across the horizon. 

If one further assumes that the Cja and Kg components of the extrinsic geometry are time-invariant {C^uja = and 
CgKf^ = 0) we have a weakly isolated horizon. Then, by (|19p (fjZbcde^'^t'^n'^ — and so gathering many of the previous 
results together, implies that rfa^f = and so the surface gravity is constant over the horizon. This is the zeroth 
law of black hole mechanics. 

There is also a phase space version of the first law of black hole mechanics which examines how physical quantities 
change across the phase space of isolated horizons [l5^. It follows from a careful Hamiltonian analysis of isolated 
horizons however it is beyond the scope of this paper. Here we will instead focus on dynamic versions of the first and 
second law, starting with event horizons. 



D. The first and second law for event horizon mechanics 



Next consider dynamic event horizons - null n-tubcs with ^ which satisfy appropriate teleological boundary 
conditions. Then the Raychaudhuri equation (|23p multiplied by the area element can be rewritten as 



Kg 



n-2 
n- 1 



(26) 



Forms of the first law of black hole mechanics may be derived from this in two distinct ways. First, the teleological 
first law. On integrating ()26p over the horizon between the slices v — vi and v = V2, it becomes 

S2 



V2 

dX 



X [ Ki + 



n-2 



n — 1 



where a is the area of the horizon. In the case where the event horizon transitions between (near) equilibrium states 
on which ct « and does so slowly (so that is small relative to nt) this gives the integrated first law 



( d\ j d^x (k^Av^) «^ dAy" ^/^(llCT^'^^f + Gahr^'') 



(28) 



As we saw in section III A( contrary to initial intution, influxes of matter or gravitational radiation do not drive 
event horizons expansions, but instead curtail existing expansions (this can also be seen in the original Raychaudhuri 
equation (j23p where an increasing matter or shear flux clearly results in a decrease in Cp9(^i^). In this version of the 
first law, the integration erases this feature and leaves behind what appears to be a more intuitive first law (first 
derived in Q|). 

That said, there is a regime where a more intuitive version of the first law holds. This is the second version of the 
first law which we will refer to as the causal first law since in this case fiuxes appear to drive the expansion. If 



then 



reduces to: 



(29) 



(30) 



and a point-by-point version of the first law (approximately) holds where the fluxes (appear to) drive the expansion. 
Note that (|^^ assumes that area rates of change are small and the second derivative of area is much smaller than 
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the first. So in cases where a derivative expansion can be made in the area element (which will be the case in section 
IIVD[) . this local version of the first law holds for event horizons. 

It is important to emphasize that both of these versions of the first law apply to event horizons that are evolving slowly. 
In the first case we assumed the expansion <^ k^i) while in the second we assumed the even more restrictive (|29|) . 
In this sense, both laws apply to horizons that are close to equilibrium. This should not come as a surprise. In regular 
thermodynamics the dE = TdS form of the first law (as opposed to a more general statement about conservation 
of energy) applies only in quasi-equilibrium situations. As such it is natural to expect the same thing for black hole 
mechanics. This will be a recurring idea in the rest of this paper. 

Finally, the Raychaudhuri equation also lies at the root of the second law. For any congruence of null geodesies we 
can alway find an affine scaling of the null vectors so that ki = 0. Then, if < and the null energy condition 
holds, implies that 

Ci6(^i) < -(1/2)0^,) (31) 

and so a congruence with < necessarily includes caustics. This is the origin of the second law of event horizon 
mechanics. Other results show that event horizons don't have caustics and therefore by the preceding argument they 
must have 0(f) > everywhere. In turn this implies that Ci^/q = \fqO(^i) > 0; event horizons have non-decreasing 
area. 



E. The first and second law of apparent horizon/FOTH mechanics 

We now turn to the physics of dynamical n-tubes which are foliated by (?(f) = surfaces. In contrast to the last 
section, we begin with the second law. 

Just as the Raychaudhuri equation (j23p was key to event horizon dynamics its generalization Eq. (jlSp is the keystone 
equation for dynamical apparent horizons. Setting (?(f) =0 we get a partial differential equation giving allowed values 
of C and if one also imposes 6'(„') < and 

CnO^i^ = Ipf - daco" - R/2 + Gabtn^ < (32) 

as motivated in section fll Al many results follow. To begin, the signature of A is strongly constrained: 



1. A is null (and so isolated) if and only if no matter or gravitational wave flux crosses the horizon: C = 
Tafc^^^ = and a'f^ = 0. 

2. A is timelike (and so decreasing in area) only if the energy conditions are violated: C < ^ Tab(."'t^ < 0. This 
can happen, for example, in the presence of Hawking radiation. 

3. Under all other circumstance, interactions with its environment cause A to be spacelike (and so expanding in 
area). In this case, the horizon is said to be dynamical and this is the second law of dynamical horizon 
mechanics. 

This set of results was first shown in Q but has since been reproved many times. In the special case of C constant 
over each S\ this can be quite easily be seen from fTS)) with the understanding that ([221) holds. More generally one 
applies a maximum principle to show these results (see, for example, 



20|) 
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Next, the first law. As for event horizon mechanics, we expect this to hold close to equilibrium and so the first 
task is to define "near" equilibrium for FOTHs. Though the intuition is fairly clear, the implementation is not so 
straightforward. A near equilibrium FOTH should be: 1) slowly expanding and 2) almost null. Neither of these is 
trivial to characterize. With respect to the first condition, the second law says that dynamical FOTHs are spacelike 
and so there is no natural flow of time along the horizon against which one can judge the rate of expansion of the 
horizon. With respect to the second condition, spacelike and null normal vectors are qualitatively different with 
squared norms that are respectively positive and zero. While one might naively say that an approach to isolation 
could be tracked by following how the norm of the normal goes to zero, this runs into difficulties. The only natural 
scaling of a spacelike vector is the unit-scaling which manifestly will not approach zero but instead become undefined if 
the normal becomes null. Thus one needs some extra structure in order to specify when a spacelike vector is "almost" 
null - either an alternative scaling or some new idea. 



Fortunately a formalism already exists that deals with these problems: slowly evolving horizons |2d. l22l l23j. This 
was developed with the twin goals of getting a better understanding of quasi-equilibrium black hole mechanics and 
developing a set of tools to track the approach (or departure) of a horizon from equilibrium. Though originally 
intended for compact horizons in regular four-dimensional spacetimes it is easily adapted to (non-compact) black 
branes in arbitrary number of dimensions. 

As noted in with O^gj — the rate of change of horizon area is related to the expansion of the ingoing null normal 
by 

^v\/^= -C6(n)^/^. (33) 

The definition of V° (|12p implies that if one relabels the foliation surfaces with A ^ A then — > n°// and 
C pC where / = dX/dX. As such, the scaling dependence of the expansion may be isolated by writing 



(34) 



with the scaling dependence restricted to ||V|| = V2C. The scaling-independent part is the expansion 6^^^ associated 
with the unit normalized evolution vector V° — V"^/||V||. Note that even though V° is not defined in the null limit. 



6*^^^ is well-defined and has the desired limiting behaviour: 6*^^^ — + as V becomes null [34| . 

Thus there is a reasonable way to characterize slow expansion. To strengthen the "almost-null" analogy one should 
also require that, to first order, the evolution of the geometry of the foliations of the horizon be characterized by 
the expansion and shear associated with l'^ (as they would be for a truly null surface). Similarly the flux should be 
approximately that of a null surface. Thus one would like to find the conditions under which 

Cvqab = {9(e) - C9(n)) qab + 2{ai'^ - Cal^l') « 2af^ , and (35) 
TabV-r'' -^„fc(r-Cn")(^'' + C7n'')«T,b«^ (36) 

(since 9(£-j vanishes). Further, one would expect these quantities to be small since on a truly isolated horizon both 
shears vanish. 

Following [20] this can all be put on a better mathematical footing by making a few assumptions about the various 
quantities appearing in Eq. (j35p . Using the notation A < B to indicate that A < koB for some fco of order one, one 
assumes that all quantities are bounded relative to some length scale J^: 

\R\ , uj'^Cja , Ma^l and \Tabtn''\ < ^ , (37) 
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and similarly derivatives of C are commensurate 

\\daC\\ < ^ and WdMW < ^ . (38) 

These are essentially assumptions that all quantities are of a "reasonable" size and the geometry not be too extreme. 
They are actually quite weak and, for example, in four-dimensions all members of the Kerr family of solutions satisfy 
the first condition if C is taken as the areal radius. 

Then, the rate of expansion is characterized by an evolution parameter e defined by: 

eV^2 Maximum C f + Tabn^n^ + 6'^„)/2) . (39) 

Note that if one wishes to allow for energy condition violating matter (so that C < is possible) then it is necessary 
to take the absolute value of the right-hand side to ensure that > 0. 

This quantity can be thought of as a generalization of 0? and is independent of the scaling of the null vectors. If 
e <C 1, Eq. can be used to show that 



ll-ii^ II < while Cllai;) II and C|0(„) | < , (40) 



and so Eq. ([351) is quantified: 



^V9a6= +(-Ce^n)qab-Cail>) . (41) 



These ideas are consolidated to get the first part of the definition of a near-equilibrium FOTH: 

Definition: Let AH be a section of a n-dimensional FOTH foliated by (non-compact) spacelike {n — l)-surfaces S\ 
so that AH — {U^S'a : Ai < A < A2}. Further let V° be an evolution vector field that generates the foliation so that 
£vA = Q;(A) for some positive function q;(A), and scale the null vectors so that V° = — Cn°. Then AH is a slowly 
expanding horizon relative to the length scale if 

(i) the evolution parameter e <C 1 where 

eV^2^ Maximum C (\\a^'''>f + Tabn'^n'' + , (42) 

(ii) \R\ and \\Cjf < 1/^2 and Tabi^'n^ < l/^^ 

(iii) (n — l)-surface derivatives of horizon fields are at most of the same order in e as the (maximum of the) original 
fields. For example, ||daC|| ^ C„iax/-^, where Cmax is largest absolute value attained by C on 5*^. 



Compared to the original definition (for three-dimensional horizons in four-dimensional spacetimes 20. l22l|). the only 
changes are in the dimension and switching from compact to non-compact horizons. A secondary change following from 
these is that the area of the (compact) two-dimensional horizon cross-sections was used in the original definition to set 
the scale For non-compact horizons this is no longer feasible and so the length scale must be set differently. For 
the black branes of this paper, the only scale is set by the cosmological constant and so we define A = — (ri— l)n/2/^2 
(the relevant scale is the AdS radius); for the coordinates that we are using is equivalent to setting ^ = 1. 
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To move from geometry to mechanics, the formahsm requires that there exist a scahng of the null vectors that satisfies 
the following (again slightly modified) conditions: 

Definition: A slowly expanding horizon is said to be slowly evolving if there exists a scaling of the null vectors such 
that C < £2 and: 

(i) ||£v'^a|| and |£v'«vl ^ e/-Sf^ and 

(ii) \CvOin^\<e/^\ 

Scaling the null vectors so that C < is motivated by a couple of considerations. First, it means that the norm 
||V|| = 2C < and so in an asymptotic approach to isolation the tangent vector becomes null in an orderly fashion 
(for example it doesn't diverge in the limit). Second with this scaling, £v "time" -derivatives will properly reflect the 
slowly evolving nature of the horizon; for example Cv-Jq = — C6'(„)^/g < e^/C and so the area expansion is also slow. 

The other conditions are motivated by the isolated horizon formalism. For a horizon in equilibrium with its sur- 
roundings, foliations can always be found so that both of these vanish and then, as we have seen, the zeroth law of 
isolated horizon mechanics directly follows. Similarly for slowly evolving, near-equilibrium horizons, these conditions 
are sufficient to enforce an approximate zeroth law: surface gravity is constant across slices to order e. The derivation 
is essentially the same as that for the true first law - one simply applies the assumptions to (j22p this time setting 
quantities to be small rather than zero. Note though that in cases of high symmetry (such as those that will be 
considered later in this paper), it often turns out that those symmetries will force daKy — exactly, independently of 
these considerations. 

Finally one can combine (fTS]) and (f20|) to show that 

KvOv = Cv6(^i) + CLv6i^n)+da{d°'C -2CGj'') (43) 
+f^M:^(V) _,_ CbVV' + 0(v)e(r)/(« - 1) , 

where t° = 1°" + Cn°'. Now for a slowly evolving horizon O^i^ and Cy6(^i^ vanish and on including the other assumptions 
and approximations this equation reduces to 

Ko/:vV^« + , (44) 

where Kq is the lowest order of the surface gravity expansion on the slice: k\i k, Kq + eni. This is the first law of 
slowly evolving FOTHs/dynamical horizons. 

F. The first and second law of alternate "horizon" mechanics 

In the previous subsections we have seen that both apparent and event horizons give rise to reasonable notions of 
black hole mechanics with the same zeroth law, similar first laws, and different second laws (in that the horizons have 
different surface areas and so different notions of entropy). Now compared to isolated horizons and stationary black 
holes, event horizons drop the — requirement while dynamical apparent horizons are no longer null. However 
both return (at least asymptotically) to become isolated horizons in the equilibrium limit. With two candidates 
in hand, it is then natural to consider other n-tubes which might act as horizons: correctly interpolating between 
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equilibrium states and obeying the laws of mechanics. Guided by the example of slowly evolving horizons we will 
continue to consider the near-equilbrium limit: essentially we will consider slowly evolving almost-horizons for which 



PI. The "horizon" n-tube should be a one-way membrane in the sense that no causal signal from a trapped surface 
should be able to cross it in the direction of infinity. 

P2. Under equilibrium conditions the n-tube should match or at least asymptote to an isolated horizon. That is, the 
usual notion of equilibrium entropy should be recovered. 

P3. Near equilibrium the n-tube should be "almost" a slowly evolving horizon (and so "almost" isolated). 

P4. Near equilibrium the surface gravity Ky should be almost constant. 

P5. Near equilibrium there should be a first law of the form (f44|) . 

P6. The area should be non-decreasing. 

P7. In the apparent and event horizon limits, the conditions should reduce the known laws for those cases. 

Guided by the experience from slowly evolving horizons, we now propose a class of n-tubes that satisfy these conditions. 
As before we will need to be careful to make sure that everything remains independent of the scaling of the null normals 
but now have the extra challenge when setting conditions that we may have ^ and C = simultaneously. Thus, 
unlike for slowly evolving horizons we can't use powers of C to remove scaling invariance. Instead we will use 0(„) 
which for standard black hole/brane solutions is non-zero and of finite size when Q^j^ = 0. 

The following structures will meet all of the guiding principles. 

Definition: Let AiJ be a section of a n-tube foliated by spacelike (n — l)-surfaces 5a so that AiJ = {UuS'a : Ai < 
A < A2}. Further let V° be an evolution vector field that generates the foliation so that L\>\ ~ a(\) for some positive 
function a(A), and scale the null vectors so that \>°- ^ ~ Cn°- . Then AiJ is a near- equilibrium n-tube (NENT) 
relative to the length scale ^ if: 

(i) the evolution parameter e <C 1 where 



(iii) (n — l)-surface derivatives of horizon fields are at most of the same order in e as the (maximum of the) original 
fields. For example, H^aC'll < Cmaxj ^ 1 where Cmax is largest absolute value attained by C on S\. 

Further the scaling of the null normals may be chosen so that 



K, 0. The guiding principles in defining slowly evolving almost-horizons will be: 




(ii) \R\ and 



2 < 1/^2 



(iv) K(v), \0(n) \ - and CnO{i) < 0, 6'(^) - Ce^„) > 



(v) \cve^,^\<ey^^ 
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(vi) ||/!v<^a|| and |£v'^v| ^ e/-Sf^ and 

(vii) |/:ve(n)l <e/^'- 

Let us consider this definition and how it meets the guiding principles. First if = on A_ff , then it reduces to 
the definition of a slowly evolving horizon. Similarly if C = it will reduce to the special case of an event horizon for 
which the point-by-point first law ([30|l holds. Further if e = then ^?(^) = and C = and we are back to an isolated 
horizon. 

We now consider where this definition differs from that of a slowly evolving horizon and the resulting implications. 
First lij contains an extra term and this has the joint purpose of ensuring that the outward null-expansion is small 
and also that e is defined even for null horizons where C = 0. Further, applying it along with conditions (ii)-(v) to 
the expression (fT8|) for Cv9^e^ one can show that 

||^W||2 < ,2/^2 j.^^^agb < ^2/^2 (45) 

in essentially the same way that these quantities were bounded for isolated and slowly evolving horizons. Note, too 
that by bounding |6'(f)^(n)| to be of the order of we retain the result that to order e the horizon will evolve as a 
d(^£-j — null surface: 

Cvqab « ait' • (46) 

If |0(^)0(„)| was of order e then this relation would be lost and a host of difficulties would follow, for example in 
deriving the first law. Similarly, to order the matter flux across the horizon will also be that for a null surface: 

r,fc VV^ « Tab^t'' . (47) 

The scaling condition 6'(„) ^ — is essentially that used in the standard coordinate descriptions of exact black 
hole solutions and also for spacelike or timelike horizons implies C ^ e^. The second law is put in by hand in the 
assumption 9(^i) — C9(n) > however note that if 0(£) > (outside a FOTH) then timelike AH with C < are 
allowed. At the same time (vi-vii) guarantee an approximate zeroth law by essentially the same arguments that we 
saw for isolated and slowly evolving horizons. Similarly, one may rerun the slowly evolving arguments by applying 
our conditions to (1^5]) to rederive the first law (|44p for these "horizons" . 

The only condition left to be considered is PI. This is clearly okay if AH is spacelike or null. Intuitively however it 
should also be okay if C < 0. As discussed in the previous paragraph, this is only possible if 6'(^) > - outside any 
FOTH. So, signals could pass out of the NENT - but only from a boundary layer outside the trapped region. This is 
similar to the situation in the membrane paradigm where the boundary is taken to be a timelike surface. 

III. GRAVITY DUAL TO BJORKEN FLOW 

A. Fluid/gravity duality and "alternate horizons" 

One of the most important recent insights within the AdS/CFT correspondence is a gravity dual formulation of 
(conformal) relativistic hydrodynamics [7|. It has been known for many years that black holes are thermal in nature. 
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In the context of gauge/gravity duality thermal states of certain strongly coupled quantum field theories have been 
understood in terms of asymptotically AdS black brane geometries. From the perspective of fluid/gravity duality 
hydrodynamics is thought of as long wave-length dynamics of non-equilibrium branes. Dissipative effects in hydro- 
dynamics imply entropy production. On the other hand, entropy on the gravity side had been linked for many years 
with event horizons and area increase theorems with the second law of thermodynamics. This leads to a unique^ 
identification of the entropy also in the context of the gravity picture of hydrodynamics in contrast with the 



phenomenological definition of the boundary entropy where some freedom remains [8l. Illf. However, as reviewed in 
the last section, there are many notions of horizon on the gravity side. These lead to distinct definitions of entropy 
which coincide in the equilibrium situation, which makes it interesting to first investigate the conditions under which 
a gravitational system is close to equilibrium. The answer can be given with the help of the slowly-evolving formalism 
generalized to the AdS case in the previous section. After identification of the near-equilibrium situations one may 
proceed to understand various notions of horizons and identify the freedom of definition of entropy in the gravity 
dual to hydrodynamics. Since the entropy is an infrared concept, its local identification on the gravity side relies on 
considerations of surfaces close to the horizon (or, in other words, in the region corresponding to the thermal scale 
on the gauge theory side). The event horizon itself is a teleological concept, and as such is sometimes hard to deal 
with. The surfaces which capture the notion of entropy are better characterized by ~ (the "alternate horizons" 
from the previous section). This uses the fact that close to equilibrium (i.e. in the case of slow evolution) apparent 
and event horizon should be close to each other. 

Understanding this in the general case might be shadowed by significant technical details. However, the fluid/gravity 
duality links solutions of relativistic conformal hydrodynamics with particular solutions of Einstein equations with 



negative cosmological constant. One such solution is boost-invariant expansion known as boost-invariant flow 24 1. 
It is complicated enough to capture many non-trivial features of the hydrodynamics and at the same time simple 
enough to enable efficient treatment of the 5-dimensional geometry in the horizon region. As such it serves as 
the a laboratory for exploring different notions of entropy and the corresponding notions in the hydrodynamics of 
4-dimensional supersymmetric Yang-Mills plasma. 



B. Bjorken flow 

Boost-invariant flow is an one-dimensional expansion of plasma with the boost symmetry along the expansion axis. If 
denotes lab frame time and x^ is the coordinate along the expansion axis, then transformation to more convenient 
coordinates - proper time r and rapidity y - takes the form 

x'^ = rcoshy, 

x^ — rsinhy. (48) 

The assumption of boost-invariance is based on Bjorken's observation that multiplicity spectra does not depend on 
rapidity in the mid-rapidity region. The most general traceless energy-momentum tensor, which is boost-invariant 



Up to horizon-boundary mapping 
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(49) 



can be expressed solely in terms of the energy density e (t) 

T,. = diag |e (r) , ~rh' (r) ~ t\ (t) ,e{r) + ^re' (r) , e (r) + ^re' (r) | 
which is a function of proper time only (in boost-invariant setup no physical quantity can depend on rapidity) 

Applicability of the hydrodynamic description means that the system is fully described by the four- velocity {Uf^u^^ = 
— 1) and temperature T{t). Boost-invariance implies that = [9r]'^ and, as expected on the basis of the form of 
the energy momentum tensor, the full dynamics is encoded in the dependence of the energy density on r. Since for a 
conformal fluid the equation of state is p = l/3e, one has 

e = eoT^ (50) 
s = \eoT^ , (51) 
where s is the usual thermodynamic entropy density. 

The general form of second order hydrodynamics equations (see [s^) is 

4 e $ 

dre = --- + -, (52) 

6 T T 

3r 3 r 2 if 

Here 3> is yy component of the shear tensor 11^'' . These equations determine the proper-time evolution of the energy 
density, which in turn determines the dependence of temperature of t. 

Substituting the form of the energy density ((50|) (together with the equation of state) into these equations leads to 
the following solution for the temperature as a function of proper-time: 



where A is a scale fixed by the initial conditions and the only arbitrary number in the construction^. In much of 
the literature (such as 

mm 

) the choice A — is made. The constants rjo, Tj^"\ A^^'' are various transport 

coefficients from the first (779) and second order viscous hydrodynamics ^. They are universal numbers related to the 
microscopic physics of the underlying quantum field theory. Their presence signals the dissipative nature of the fiow 
- the entropy production. 



C. Validity of the hydrodynamic description 

The modern view of hydrodynamics is similar to that of effective field theory. It is a phenomenological description 
of phenomena on scales much larger than those of their microscopic dynamics, constructed as a systematic expansion 



It is easy to understand the presence of A if one considers the dilatation x^^ — > ax'^. Then r — > ar and A — > a~^^^ which leads to 
e (ar) = a~*e (r) which is the correct scaling 
® More precisely these are dimensionless numbers related to the transport coefficients by the relations 



Ai = X^eoT^ (54) 



where eO is defined as previously by e = cqT'' 
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in gradients. In the context of boost-invariant flow this translates into an expansion in powers of 1/t2/3. As in the 
case of any perturbative expansion, one needs to observe the regime where the expansion can be reasonably expected 
to apply. A criterion for this is that the subleading terms in the expansion be smaller than the leading order. In 
the context of Bjorken flow this can be understood as a condition on the minimal time when the expansion can be 
trusted. For the first subleading term in (|55|) to be smaller than the leading order^ one needs r > r^m, where 

TTT ■ ~;= — ; — = a < 1 (55) 

mm 

One requires that the expansion of any physical quantity such as energy or entropy density should have this property. 

Note that including higher order terms does not extend the regime of validity of the hydrodynamic expansion, but 
rather improves the accuracy within the hydrodynamic window. Moreover, in a general boost-invariant dynamical 
situation Tmin does not coincide with thermalization time, since the non-hydrodynamic (exponential) modes do not 
have to be negligible at that time [36 1. 



D. Holographic description of Bjorken flow 

Since the energy- momentum of the gauge theory is related by AdS/CFT dictionary to the five-dimensional asymptoti- 
cally AdS metric, three independent warp factors are to be expected on the gravity side. This leads to a boost-invariant 
metric ansatz of the form 

ds^ = 2g (f , r) df dr + h (f, r) dr^ - r^A {f, r) df^ + (1 + rff e^^^'^'My^ + r^e^^^-'^Mxi (56) 

The functions g{f,r) and h{f,r) are gauge degrees of freedom and can be chosen conveniently. Setting g(f,r) to 
and h (f, r) to 1/r^ leads to the Fefferman- Graham-like coordinates, which are particularly useful in obtaining the 
boundary energy-momentum tensor. These coordinates suffer however from an important drawback: they break down 
at the locus where A (f , r) goes to zero, which is where one would naively expect the location of the horizon. In order 
to keep a full control over the geometry in this region it is useful to adopt so called ingoing Eddington - Finkelstein 
coordinates which give the following metric ansatz |26| 

ds^ = 2dfdr - r^A (f , r) df^ (1 + rff e^^^""'^dy'^ + r"^ e"^^ '''^ dxi (57) 

This form of the metric still has a residual gauge freedom r r + R{T) In the following this freedom will usually 
not be fixed, since it provides a useful cross-check on the correctness of the calculations - physical quantities in the 
boundary theory should not depend on the choice of gauge in the bulk. 

The metric (j57p describes the gravity dual to boost-invariant plasma for any r. It turns out however, that in the 
regime of large proper time, the boundary energy density e (r) satisfies the equations of hydrodynamics. This highly 
nontrivial observation was first made in [2'| and then developed in using the language of gravity dual. 

It is particularly easy to understand r^^/"^ damping of the subleading pieces in the expression of the temperature - 
they correspond to the gradient expansion. The gradient expansion has to have its counterpart in the gravity dual to 



The temperature has been chosen here because it enters the definition of the gradient expansion. 
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boost-invariant flow in the regime of large proper time. The boundary gradient expansion leads to the introduction 
of a scaling variable and the large proper time expansion in the scaling limit (keeping v — r ■ t^I'^ fixed while taking 
the limit f — s- cx)) on the gravity side 



C (f , .) = CO (. f V3) + J_,^ ,1/3) + ,1/3) + . . . . (53) 



The meaning of the scaling limit is that the radial distance in AdS is measured in units of inverse temperature 
{r/T ^ r ■ T^^/"^ since T ^ t^^/'^ in the leading order). The large proper time expansion corresponds to the gradient 
expansion in the hydrodynamics. The metric with Aq (rf^^'^) only will be called the zeroth order metric and is dual 
to perfect fluid hydrodynamics on the boundary side. The first order gravity solution - Ai (rf^^^^, bi (rf^/'^) and 
ci (r f ^Z"^) - mimics boundary viscous effects on the gravity side and is the leading order relevant for the gravitational 
entropy production, whereas the second order metric corresponds to relaxation time in the boundary theory. The 
formulae providing the form of the warp factors at zeroth, first and second orders are given below, whereas third order 
formulae can be found in Issl . 



At order zero one has 



4a4 



TT^A 



Ao{v) = 1 
bo{v) = 0, 

co{v) = 0. (59) 



which is equivalent to the Janik-Peschanski solution For simplicity, this solution assumes a choice of gauge to 
bring out the similarity to the static AdS-Schwarzschild black brane (see section HVEp . 



At first order one finds 

Ai (v) 



27r3(« + 7rA)A3 2 {v^ + tt^A^) 6^ 



biiv) = - 



2arctan(^) 41og(i;) , log + 21og(^ + l) 



SttA SttA 37rA 37rA 

41og(^) 2Ji 4 1 
SttA 3t; 3v 3A' 



ciiv) = -l6i(«) + ^. (60) 
2 V 

The second order formulas can be found in appendix |B] 

The zeroth order solution looks like the boosted and dilated black brane. This is in line with the argument of which 
states that the gravity dual to the general solution of hydrodynamics can be obtained by boosting and dilating black 
brane and having this as leading order solution solving Einstein equation perturbatively in the gradient expansion. 
Since boost-invariant flow is a particular solution of the equations of hydrodynamics, this is not a coincidence. 
The interesting, yet not surprising given its interpretation as one-dimensional expansion, feature of boost-invariant 
hydrodynamics is that the temperature does not settle down to the constant non-zero value, which gives rise to the 
question as to whether this solution is covered by an event horizon. In the next sections it will be demonstrated that 
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this indeed the case, which means that the gravity dual to boost-invariant flow is a well defined geometry from the 
cosmic censorship hypothesis point of view. 

The gravity solution displayed above determines, via holographic renormalization, the boundary energy momentum 
tensor [s^, which is of the form (j49|) . Explicitly, one finds 

, , 3 2 2 2 1 , l + 21og(2) 1 ^ 



8 " T^l^ \ SttA t2/3 187r2A2 r^/a 

-3 + 27r2 + 241og(2) - 241og^(2) 

4867r3A3T2 ^' 



(61) 



Furthermore, the proper-time dependence is exactly what is required by hydrodynamics. The temperature (determined 
by (HI]) and ^) reads 



1/3 \ 67rA r2/3 367r2A2 t^/^ 
-21 + 27r2 + 511og(2)-241og2(2) 1 



19447r3A3 t2 ' ^^^^ 

Comparison with the static situation determines the coefficient eo in ([50)1 to be 

eo = ^NItt^ . (63) 

Matching the energy density with (jSOp and (|53p and using the linearized hydrodynamics to obtain the relaxation time 
from sound wave dispersion relation determines the transport coefficients of N=4 Yang-Mills plasma: 

^ = (64) 
2 - log 2 

Ai = ^ (65) 



These expressions are based on terms up to second order. The formulae (|6ip and ([621) include third order terms, but 
since the tensorial structure of the third order hydrodynamics has not been yet investigated ^'^ , the relevant transport 
coefficients cannot be calculated. 



IV. HORIZONS IN THE BOOST-INVARIANT SPACETIME 



A. Preliminaries 



With the tools from the last section in hand one can now turn to the identification and study of various notions 
of horizons in the spacetime defined by the bulk metric ([57]) . Since the main goal is to study possible notions of 
entropy in the dual field theory, the hypersurfaces of interest are those which satisfy the symmetries of the boundary 



It seems that there is no practical need to do so - the inclusion of second order hydrodynamic terms in RHIC simulations changes the 
results by just couple of percent [loj. The third order terms will give even more suppressed contribution. However there are hints that 
in some cases resummation of hydrodynamic might be needed |4l| . 
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dynamics under consideration (Bjorken flow). This singles out spacelike three-surfaces of constant r and r. Such 
surfaces possess outward and inward pointing null normals 



+ ^r^A{T, r) 



d_ 

dr 



n 



d_ 

dr 



(66) 

-1. The rescaling freedom of ^ — > /(t)£ and 



As usual, these vectors have been cross-normalized so that I ■ n 
n — > n/ /(t) remains but the specific scaling used here is chosen to be consistent with the flow of time at asymptotic 
infinity. This will be discussed in more detail below. 

The hypersurfaces of interest will all lie within the r = constant hypersurfaces. Given that coefficients of the metric 
are given by the series (I58p , solutions can be sought in the form 



1 



iro 



1 



1 



The coefficients rt appearing here will be determined by the conditions imposed, and it turns out that the solutions 
are unique. Not only is the event horizon uniquely defined, but in this case demanding that the FOTH shares the 
symmetries of the spacetime means that we can also select a unique FOTH. In this section the focus will be on 
the FOTH and the event horizon, while in the following section more general surfaces of the form (j67p will play an 
essential role. 



rg^...). 



(67) 



B. The boost-invariant FOTH 



To find marginally trapped three-surfaces within the r = constant hypersurfaces, one needs to solve the equation 
6'(^) = 0. Evaluating 6'(£) on a hypersurface of the form ([67)) yields 

/ • 

Solving fl^^-j — (to third order) shows that there is a unique apparent horizon on the r = constant slices at: 

, , 1 ( 1 f Si 1\ , f fAS2 , 1 1 log(2) 

rAH{T) = -jT^T <^ vrA ' 



(69) 



7-1/3 I 7-2/3 \^ 3 2 J Ar4/3 ^ 3 97r 24 IStt 
f / C A63 25 1 1 71og2(2) 



A2t2 Vl87r2 3 4327r 8f7r2 7776 f627r2 
log(7rA) 2 log(7rA) 25 log(2) 



IStt 277r2 1627r2 

For the three-surface defined by r = rAH{T), the inward expansion is 

^ 1/3 f 3 1 -l + 21og(2) 1 / 7 log" (2) log(2) 

^ \ ttA 27r2A2 t2/3 127r3A3 r^/a A^ t2 547r4 2471^ 

251og(2) 1 C 35 Y 

547r4 2592^ ^ 6^ 216^ j 



(70) 



where C is Catalan's constant. As expected, (|70)) is independent of the constants 6k which reflect the gauge dependence 
of (|57)) . It is easy to see numerically that 

0.95 0.051 0.0010 0.00039 

V) - J- - ^2^2/3 ^ AVV3 ^ A4t2 ^ ^ 
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This is clearly negative and will stay negative in a neighbourhood of tah- Further it is clear from the expression for 
the outward expansion that for r « rAH, r > tah ^ (^(e) > and r < tah ^ (^(t) < 0. That is, there are fully 
trapped surfaces "just-inside" r = rAH and so this marginally trapped surface bounds a fully trapped region and so 
can be identified as a black brane apparent horizon or FOTH. 

The remaining geometric quantities discussed in section |TT] can now be calculated. First, requiring that the evolution 
vector V° be tangent to the horizon, implies that in the large t expansion the expansion parameter C has the form 

C = C-i + ^Co + ^C, + ^C, + -L^C, + 0^-^^ , (72) 

for some set of coefRcients C_i, . . . C3. In this case it is straightforward to see that C^i — identically in consequence 
of the structure of the large r expansion (even without using the explicit form of the solution). The coefficients Cq 
(perfect fluid) and Ci (viscosity) also turn out to vanish, so the leading contribution appears at order All in all, 
on the r = rAH{T) FOTH one finds 

C=(iU-^(i|ei l),o(^V ,73, 



9 J r2 Ar8/3 \ 9tt 54/ ^^n/a 

This is clearly greater than zero and so the horizon is dynamical: that is spacelike and expanding. This can then be 
cross-checked in two ways. First one can directly calculate the volume element on the three-surfaces: 

e = \/^dy A dxi A dx2 (74) 

where, up to third order 



1 1 o .0 A /I 



4f4-T^-:^ + ^-i-log(2)-^^^ 



^216 144 2592 216 24 °' ' 216 
Then, to lowest order, the rate of expansion is 



1 _ J_ (J_ 1 f log(2) _ J 



dT t5/3 \ SttA A2 t2/3 1^ 37r2 IStt IStt 
which is clearly positive. Alternatively 

1 d^ 1 



(76) 



^£vV'? = -C0(„), (77) 



y/q dr 

and substituting in the appropriate values obtains the same result. 
Other quantities arc the squares of shears in the two null directions: 

1 /351og2(2) log(2) 141og(2) 1 1 1 



A2ri"/3 V 547r2 67r 277r2 277r 97r2 648 

An)^a, ^ J 3_ _J_ ( 31og(2) l\ , 

^afc ^(n) ^2/3 8^4A4 + ^5^4/3 1^ 16^4 + 2^5/ ^ ' 

1 / 351og^(2) 31og(2) 51og(2) 1 1 25 

A6r2 I 967r6 327r5 dn^ 11527r4 Stt^ A8tt^ 
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and the accelerations/inafRnities in the two null directions: 

= --^'V.r = ^ {2.A - + ^ (i. - MS I ,80, 

1 f7\og^{2) log(2) 71og(2) 1 C 1 1 



A2^4/3 \^ 8l7r2 367r 8l7r2 547r 97r2 1Q2tt^ 3888, 

K(„) = -£an''Vbn'' = (81) 

Since vanishes we have 

«(v) = «(^) - C" = • (82) 
Finally, the connection on the normal bundle is 

LVa = -g^n^Vb^ = , (83) 
and the (three-dimensional) Ricci scalar of the r = constant and v — vah three-surfaces is 

R = 0. (84) 



C. Boost-invariant flow and slow evolution 



Having described the geometry of the FOTH in the boost-invariant spacetime, the next step is to try to understand 
its physics. Given that the metric (j56p is a perturbation of a boosted black brane and that a boosted brane is simply 
a coordinate transformation of a static black brane spacetime, one would intuitively expact that the horizon should 
be in a near-equilibrium state. Thus the physics of these horizons should be quasi-equilibrium physics. Quantifying 
this intuition and then understanding its implications will be the subject of this section. 

One can easily check that these conditions hold for the FOTH in the boost invariant black brane spacetimes. First, 

^^-2.^A^7^(l + o(-i,)) (85) 

and the evolution parameter for the horizon is of the same order as the expansion parameter for the metric. Next (ii) 
and (iii) are trivially seen to hold: there is no dependence of the geometry or scaling of the null normals on {w, x, y) 
and all of these quantities vanish. Thus the horizon is slowly expanding to order l/r^/"^. Though it is a consequence 
of these conditions, in this case we can also explicitly check that the horizon is almost null in the sense of Eq. (|35p : 
the evolution of the three-slices is characterized by the expansion and shear associated with just as they would be 
for a truly null surface. 

One can now check these conditions (and their consequences) for the null scaling ((66|) . Though it is clear that the 
conditions will hold for a variety of scalings, we make this particular choice so that "time" evolution on the horizon 
is consistent with that on the boundary (both are proportional to r) and, as a retroactive justification, the surface 
gravity matches the temperature on CFT side. Explicitly checking the conditions, the first part of (i) holds trivially 
since (Dq = by equation (l83|) while 

Kv = -VribVaf^ = (86) 
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by equation ((8T|) and so from (|80p the second part holds as well: 



l^v^^v ~ -^"^^ =^ \^vKv\ < . (87) 
Finally for (ii) a similarly trival calculation from ([70]) shows that 

^vO(n} ~ =^ \^vO{n)\ < e- (88) 

Thus, the horizon is slowly evolving to order e ~ r^^/'^. Then, in addition to it being geometrically "almost" null it is 
also mechanically close to equilibrium. As already noted, Einstein's equations are sufRcient to imply that the zeroth 
law of black hole mechanics (almost) holds: the surface gravity k\> is approximately constant on each slice. Here, the 
symmetry of the black brane means that a stronger result is valid ~ from ([50]) one can see that the surface gravity is 
constant on each slice (though it evolves slowly up the horizon). The dynamical first law also automatically holds, 
but again in this case one can also check it explicitly. To order r^^ one has: 

^^v/:vv^«y^||a(^)||'«|^'A3l, (89) 

O 7" 

and so, as expected, the expansion is driven by the null shears (which would normally suggest incoming gravitational 
radiation). The dependence on t in this expression matches what is expected in leading order on the basis of 
thermodynamics of Yang-Mills plasma. To see this, note that the volume in the boost invariant setting depends 
linearly on r, so that the entropy of the plasma can be written as 5* = stVq where Vq is a reference volume and s is 
the entropy density ([5T|) . Using the results of the previous section one then finds 

which is consistent with 



D. Event horizon 



It is usually said that the event horizon is somewhat inconvenient to work with, since determining it requires knowing 
the entire future evolution of the spacetime under consideration. This is indeed an onerous requirement in the typical 
situation of computing the evolution of spacetime geometry starting from some initial data. The setting explored in 
this paper is in a sense complementary: the spacetime geometry is constructed order by order in a large proper-time 
expansion (or gradient expansion) starting in the far future at zeroth order. This circumstance makes it possible to 
determine the location of the event horizon in the late time regime. 

The method of finding the event horizon for boost-invariant flow closely resembles the one presented in Q for the 
gravity duals to fluid dynamics Q. The crucial assumption there was that the metric relaxes to (uniformly boosted) 
AdS Schwarzschild, where the position of the event horizon is well known. The event horizon for the metrics there was 
defined as a unique null surface which asymptotically coincides with the event horizon of the static AdS Schwarzschild 
dual to the uniform flow at constant temperature. Despite the fact that this is not the case for boost-invariant flow, 
it is still possible to find a unique null surface which can be interpreted as an event horizon. The key observation 



Obviously leading order means first order in the gradient expansion, since entropy is preserved in the perfect fluid case. 
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is that the event horizon should coincide with the FOTH in the large-proper time regime and within the scaling 
limit^2 jts radial position in AdS should depend on proper time only, which reflects the boost-invariance (no rapidity 
dependence) together with translational and rotational symmetry in the perpendicular directions (no x± dependence) . 
If r is the radial direction in AdS space, t the proper time and teh (t) expresses the time evolution of the horizon, 
then the equation defining the sought co-dimention one surface in AdS takes the form 

r - TEH (r) = 0. (91) 

The covector normal to the surface is dr — (t) dr and requiring it is null gives the equation for teh (t) 

AiT,rEH)-rlH-2r'EH = 0, (92) 

where for clarity the dependence of teh on t is ommited. This equation can be solved perturbatively in the scaling 
limit. 



Using the solution described in [26] (valid up to third order in the late proper-time expansion) one finds 

1 , fS2 , 1 1 log(2)\ 1 , 
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(93) 



The dependence on arbitrary constants of integration reflects the gauge freedom - the position of the horizon is a 
gauge-dependent quantity {r ^ r + f (t) for arbitrary / (t)). Comparing with ([M]) it turns out that the FOTH 
coincides with the event horizon in the leading and first subleading orders. This is in agreement with the observation 
that the constant C is non-zero only in the second and higher orders in expansion. The second orders differ and 
the FOTH becomes spatial. 

It is be also interesting to compare the position of the apparent and event horizon with the naive horizon defined as 
the hypersurface on which A(t, r) vanishes 



■n^'™ - ^1/3 Y'^- ^2/3 \^3 3 , 

Such a hypersurface coincides with the event and apparent horizon in the leading order, but differs in the first 
subleading order - it is situated between the event horizon and the boundary. This means, for example, that truncating 
numerical simulations at the point where A{t, r) vanishes is simply wrong. 

Finally, observe that in the naive limit t oo the boost-invariant metric relaxes to the empty AdS5 metric instead 
of the static AdS-Schwarzschild solution. However, this is not so strange from the dual CFT point of view, where 
the fluid is expanding to infinity and its energy density becomes smaller and smaller. It means that the boundary 
system does not permanently thermalize to non-zero temperature. The interesting feature of the boost-invariant flow 
is an apparent thermalization, which expresses itself as an applicability of the equations of hydrodynamics in the late 
stages of the evolution. 



The scaling limit insures that the higher derivative corrections to the geometry are small (MH thanks Toby Wiseman for discussion on 
this point) 



30 



E. Revisiting the scaling limit 

The symmetries of boost-invariant flow make it possible to seek the location of the event horizon considering only the 
variables r and t. It is possible then to focus only on the dr, dr part of the full metric, which at leading order takes 
the form 

ds^ = 2drdr - |l - dr' + . . . (95) 

The scaling limit discussed at length in section Hill involved introducing the scaling variable v — rr^^^ which is kept 
fixed as r — > cx). This motivates the following change of variables 

3 

which leads to 



(96) 



ds^ = 2dudv - ^2 |l - du'^ + O (^i^ . . . (97) 

This shows that this part of the metric looks precisely the same as the corresponding part of the static black brane 
metric, with v denoting the radial coordinate and u Eddington-Finkclstein ingoing time coordinate. This means that 
at leading order in the late-time expansion the problem of determining radial geodesics^^ in the asymptotic boost- 
invariant geometry is the same as in the static case. It is then not surprising that the naive position of the horizon 
coincides asymptotically with the actual event horizon. 

Note however that these considerations do not imply that the asymptotic geometry is static. Clearly, the remaining 
terms in the metric are time-dependent even after this coordinate transformation, even though the area of the event 
horizon remains constant in this regime. 



V. PHENOMENOLOGICAL NOTIONS OF ENTROPY 



A. Introduction 



It is believed that equilibrium states of black holes are thermodynamic in nature. Their entropy is associated with the 
area of spacelike slices of the event horizon in an unambiguous way and the second law of thermodynamics is linked 
with Raychaudhuri's equation (|23p . The property that the area of the event horizon is non-decreasing continues to 
hold in a generic dynamical setting. This prompts the question whether there is a sensible notion of entropy valid in 
such a non-equilibrium situation. However, as discussed in section (|IIF|) . there are hypersurfaces of non-decreasing 
area other than the event horizon (which coincide with it in the static case). The notion of entropy thus becomes 
less clear in these cases, as frequently discussed in the literature. The AdS/CFT correspondence makes it possible to 



As stressed previously, this is all that is needed to determine the location of the event horizon. 
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view this problem from the gauge theory perspective. For example, the teleological nature of the event horizon leads 
to acausal behaviour of gauge theory entropy associated with it 12 1 . 

The bulk description is under control precisely when the field theory is strongly coupled, which in itself makes it 
hard to analyse directly. However it is well understood how to formulate a hydrodynamic description of the system 
in the appropriate regime. This description depends only on symmetries and the idea of the gradient expansion, as 
explained in 0, [s^ for the conformal case and in [ll]| in general. 



B. Entropy from second order hydrodynamics 

The requirement that entropy should be non-decreasing during hydrodynamic evolution can be expressed in a covariant 



way in terms of an entropy current whose divergence is non- negative 42] ■ While the energy momentum tensor is a 
canonically defined operator, the entropy current is a derived notion. In the spirit of hydrodynamics (or effective field 
theory) it is also constructed in a gradient expansion as the sum of all possible terms at a given order. The dynamical 
equations of hydrodynamics are the conservation equations for the expectation value of the energy-momentum tensor. 
Thus, the coefficients appearing in the gradient expansion of the expectation value of the energy-momentum tensor 
(the transport coefficients) are the physical parameters of this phenomenological theory, since they figure directly 
in the evolution equations. They describe physical properties of the underlying quantum field theory. In contrast, 
the coefficients which appear in the phenomenological expression for the entropy current are constrained only by 
the requirement that its divergence be non-negative. These parameters are logically independent of the transport 
coefficients. At the present level of understanding they reflect a real ambiguity in the phenomenological notion 
of entropy current in hydrodynamics (as explained in the following subsection). This ambiguity is however of no 
consequence when entropy differences between equilibrium states are considered. 



In the case of conformal fluids the most general form of the entropy current was recently constructed [8|, |ll| up 



to second order in gradients. The crucial symmetry requirement was Weyl covariance (see 43| for useful explicitly 



Weyl-covariant formulation of hydrodynamics). Using the notation of [ll^ (which descends from [8|) one has 

+Bi i -Vf + —Sij + B2 {V^ - u'^52) . (98) 

Here s denotes the thermodynamic entropy density (|?T|) . and 5i_2,3 are the 3 possible conformal (Weyl-covariant) 
scalars 

Si = a^.a''- , 52 = n^M'''' , 
53 = c2 V;^ V'[ In s + In sV'^ In s - ^Uc,upR°''^ ~~\R+l(y-uf, (99) 



,14 



In the hydrodynamic formulas A**" = is projector to the fluid's local rest frame, whereas = A*"^ V^. Moreover fluid shear 

tensor (responsible for dissipation in the first order conformal hydrodynamics) reads ct'"' = A''" A"'' (Va«^ -I- V^m^) — ^ij- A^" Vat*" 
and vorticity (nonzero for rotating fluid) = A''°'A'''^ (Vq«^ — V^Mq). 
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and Vi^2 are 2 possible conformal vectors^^ 



where a is the hydrodynamic shear tensor and f2 is the vorticity, and R is the boundary Ricci tensor. 
The entropy current depends (|98p on 5 constants Ai 2,^ and -B1.2 and its divergence reads 
-^An-.^ = iv>V(-2A3 + i3i) + iv^a^''V^hrs(-2A3 + Si) 
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(100) 



(101) 
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If the shear tensor is non-vanishing^^ the positivity of the shear viscosity r/ should guarantee (see however llj and 
the next footnote) that divergence of the entropy current is non-negative: higher order terms cannot change this 
conclusion as long as the gradient expansion is valid. However, as noted in 8], it is perfectly reasonable for a^^ to 
locally vanish (requiring this imposes just 5 conditions for derivatives of the four-velocity) in which case the higher 
order terms will dominate the entropy production. Positivity of poip thus requires 



= 2^3 



(102) 



Since the shear tensor a^^ is multiplying the whole square bracket in (|10ip . in the case when it is the whole 
contribution from first two orders is absent. At this level there is a real 4-parameter ambiguity in the hydrodynamic 
construction of the entropy current^^. 



C. The case of Bjorken flow 



From the perspective of the AdS/CFT correspondence it is natural to ask whether the ambiguities appearing in the 
construction of the hydrodynamic entropy current match on both sides of the duality. In a general situation this 
might be involved, but one may try to gain some insight into this question by considering a particular solution. The 
Bjorken flow provides a simple, highly symmetric, yet non-trivial example. 

The current l|98p evaluated on the boost-invariant solution given by the velocity = \dr\^ and temperature T(r) 
((53)) takes the form 

= ~au^ (103) 



Note the absence of parity breaking terms present in [^. For discussion of hydrodynamics with parity breaking terms see [44l . 
Note that cr^iyO"^" as a trace of the square of matrix cannot be negative and vanishes if and only if u*^^ = 0. 
^"^ Considerations of the case with an arbitrary small tr'"' in suggest that further constraints on the entropy current may be imposed. 
In particular, the only freedom left after such consideration is in the parameter Ai . It appears that these arguments rest on competition 
between terms of different orders in the gradient expansion. This paper cannot shed any light on this, since in the boost invariant case 
the shear is not close to vanishing. 
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Sir) = s(r(r)) ^1 + 2 } (104) 



with 

37r2r(T)2r2 

In general the entropy current does not have to be proportional to the flow velocity beyond leading order (perfect 
fluid), but in the special case of boost-invariant flow non- leading order effects are captured by the single scalar function 
s(r). This function involves 3 of the 5 constants appearing in the general phenomenological construction. Changing 
the value of Ai has been identified with the freedom in choosing the horizon to boundary map [8]. The ambiguity 
parametrized by A3 was not interpreted in [sl . One would like to interpret this freedom in terms of allowed definitions 
of "horizon" on gravity side. Note that the example of Bjorken flow, while rather special, is still rich enough to 
partially capture this ambiguity. 

In quantitative terms this ambiguity can be estimated as follows. In order for the hydrodynamic expansion of p04p 
to be valid the magnitude of |Ai + A3I should bounded so that the leading term dominates for times larger than Tmin 
defined by (|55p . Expanding (|104[) up to second order one has 

~. . A^r 11 (8(Ai+A3)+log(2)) 1 1 

^ V 1' - 2^757^ + 12^5X5 ^\ (105) 

Demanding that the second order contribution be smaller than the first order correction by a factor of aP at r = Tmin 
leads to the bound 

i(-/3 - log(2)) <A,+A,<^{P- log(2)) (106) 

where /3 is at most of order 1/a. This provides a rough estimate of the allowed inderminacy in the phenomenological 
notion of non-equilibrium entropy as defined by (1104^ : 

One would like to understand this quantitively in terms of the freedom of defining entropy on the gravity side. 



D. Entropy from gravity 

As reviewed in section|lTl in a dynamical setting it is no longer clear if there is an appropriate geometrical notion which 
can be used for the definition of entropy. For example both apparent and event horizons appear to give rise to notions 
of entropy, which satisfy the second law of thermodynamics and coincide in equilibrium. This provides motivatation 
to look more generally at the dynamics of hypersurfaces whose area is non-decreasing. The starting point should be 
equation (jl3|) . This formula determines the rate of change of the area element of a general 3-hypersurface in the terms 
of the expansions 0(f) and the expansion parameter C . In the boost-invariant case the 3-hypersurfaces consistent 
with the boundary symmetry have the form (j67p . Their area is given by 



where 



ri=ri- -di 

h = r2 + ^2 (109) 



34 



The notion of entropy defined by such hypersurfaces is 



Ar2 

S=-^A (110) 
zvr 



(since in the units used here (AdS radius set to 1) = 27r ^N^)- 

Equation (|13p , which expresses the change of area of the hypersurface sections becomes 

6f 1 + 2 1 18f ? + 12f 1 - 27r(18f2 + 1) + 6 - 5 log(2) 1 



SttA f5/3 97r2A2 f7/3 



(111) 



In this formula ro has been set to tt to match the thermodynamic entropy when aU gradient corrections are discarded. 
For the leading term (at order 1/f^/'^) to be non-negative one gets the bound fi < —1/3, and then requiring that the 
following term be smaller gives an allowed range for r2- At this level of analysis this is all one gets; ri is not fixed^^. 
Both notions of entropy (defined via the event or apparent horizon) provide unique ri, which lies in the allowed range. 

The entropy density^^ obtained from the third order expression for the event horizon reads 



2.3 [,11, 6 + 7r + 61og(2) 1 



■^EH (r) = -A^.^vr^A^ 1 



2 " i 27rA f2/3 247r2A2 f4/3 

420 + 907r + 7r2 + 372 log(2) + IOStt log(2) + 420 log^ (2) 
25927r3A3f2 

In the dynamical case the result for the entropy density reads 



(112) 



2 " T \ f 2/3 2nA A2f4/3 \^ 47r2 24n 12n^ 
1 /351og2(2) log(2) 51og(2) 1 1 



A3f2 V 2167r3 247r2 2167r3 25927r 1447r2 2167r3 



(113) 



Numerically one finds 



For the event horizon 



l,r2 2,3lf, 1 0.16 1 0.039 1 0.00651 



1..2 2a31 0.16 0.056 0.0181 

The key observation is that the event and apparent horizons coincide at the leading and first subleading orders, 
which is a hint that also in the case of a general surface there should be no ambiguity until the second subleading 
order. If one is to identify the entropy defined here with the field theory observable, as required by the AdS/CFT 
correspondence, then it should be Weyl covariant in the boundary sense. To do this explicitly one would need to solve 
the Einstein equations with the boundary metric given by 

dsl = e-^-'f^) {-dr^ + T^dy^ + dxl } (116) 



Note in particular that there are surfaces outside the event horizon which are acceptable from this point of view. 

Entropy density is understood as entropy per unit volume, which in the proper-time - rapidity coordinates involves a factor of r. 
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where w(t) is a conformal factor having the form of an expansion in powers of The entropy computed this 

way would be Weyl covariant (i.e. proportional to the appropriate power of the conformal factor) only for fi = —1/2., 
i.e. the value assumed by fi in the case of the event or apparent horizon (which coincide at this order). The quick 
way to get this answer is to write the entropy ()110|) in terms of temperature and velocity, whose transformation rules 
under Weyl rescalings are known. This procedure parallels the field theory analysis reviewed earlier. The first step 
is factor out the thermodynamic entropy which sets the Weyl transformation property of the entropy density. This 
leads to 

, 1^^2_2^.,^3 , 6f i + 3 , SOff + 42^1 + 367rf 2 + 2^ + 9 + 4 log(2) ] 

Since there are no Weyl-covariant scalars nor vectors at first order in derivatives, the only way that this formula can 
be Weyl covariant is if the first subleading term vanishes, which determines ri = ^1/2. When this result is used in 
eq. (|llip . one finds 



^v\/q = \/l{0(i>) - C6'(„)) 

1 1 9 - 47r(18f2 + 1) - 101og(2) 1 



37rAf5/3 187r2A2 f7/3 



(118) 



The leading contribution ensures positive entropy production due to the shear viscosity so there are no further con- 
straints on f2 . The appearance of possible freedom in choosing f2 can be understood following again the hydrodynamic 
argument presented in subsection (|V Bp . In the first place, note that the formula (|118p is evaluated using the bulk 
Eddington-Finkelstein proper-time f, which raises the question about its relation with boundary proper-time coordi- 
nate T. Such a mapping freedom has been addressed using Weyl-covariant language in Q and amounts to the trivial 
mapping in the first two orders of the gradient expansion, with an ambiguity showing up at the second order. In the 
case of boost-invariant flow the most general mapping (up to second order) takes the form 

- — -d+A^) (119) 

where 5Ai is a constant parameter multiplying Weyl-covariant scalar Si. For such mappings to make sense within 
the context of the gradient expansion this parameter must be suitably bounded as explained earlier. This shifts the 
second order contribution in (|114p and (|115p . so it must be of order in the sense of ([55]) . which leads to 

SAi . ±^ (120) 

with 7 of order 1. This mapping freedom accounts only for a part of the ambiguity in the phenomenological construc- 
tion of the hydrodynamic entropy current. The rest of the ambiguity can be understood following the Weyl analysis 
of the properties of gravitational entropy. At second order in gradients there are two Weyl covariant scalars {Si and 
S^) and one vector (Vi) which are non-trivial when evaluated on the boost-invariant solution, with the vector being 
proportional to velocity in that particular case. Since mapping freedom is identified with the Si contribution^" , it is 
clear the f2 must come from the relevant combination of S^ and Uf^V^. Evaluating (|117p with fi = —1/2 gives 

lAr2 2^^ n3 f-, , -3 + 367rf2 + 27r + 41og(2) 



-«,Vr(rni+ (121) 



It is also reasonable to stipulate that Ai is partly linked with the mapping freedom and partly with a freedom in choosing between 
different notions of horizons in the bulk. 
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Comparing this with (|104p one finds 

Although this is all one can get from the analysis of the gravity dual to the boost-invariant flow^^ it is reassuring 
that at least in this case the gravity picture is capable of capturing the ambiguities of the boundary phenomenological 
construction. 



E. A phenomenological definition of black brane entropy 



The freedom in the definition of the hydrodynamic entropy current on the gravity side follows not from the various 
possible notions of horizon, but rather from adopting the phenomenological construction in the bulk, which is analogous 
to the boundary one. The surfaces considered are not horizons in any of the usual senses, but they do have the property 
that their area increases, and that is all that is required if one adopts the slowly-evolving geometry approach. In order 
to understand this in greater generality one would of course need to go beyond the Bjorken flow example and consider 
the equation (|13[) evaluated on the gravity dual to the general hydrodynamics. This is certainly possible employing 
the Weyl-covariant formulation in the bulk. It would be very interesting see this in detail. 



VI. CONCLUSIONS AND OUTLOOK 



The main goal of this article was to explore the relationship between the notions of entropy on both sides of AdS/CFT 
duality. This lead to a phenomenological definition of black brane entropy on the gravity side, which was inspired by 
the corresponding construction in hydrodynamics. In the case of Bjorken flow the freedom inherent in this definition 
accounts for the entire ambiguity appearing in the hydrodynamic entropy current in this case. This leads to an 
understanding why the event horizon coincides in the leading and first non-leading orders of the gradient expansion 
with the unique apparent horizon compatible with the boundary flow. The origin of this circumstance is the Weyl 
invariance of the boundary theory. It would be very interesting to understand the form of general hydrodynamic 
entropy current from the bulk perspective and carefully understand the sources of possible freedom in such dcflnition. 

It is natural to ask what is the physical relevance of the potential ambiguity in the definition of entropy current. In 
the case of local entropy production such an ambiguity signals lack of physical meaning. This however should not be 
disturbing, because the thermodynamic notion of entropy makes sense only in equilibrium. Since one expects that 
systems described by hydrodynamics equilibrate due to dissipative effects, the total entropy can be calculated in the 
late stages of evolution and is given by the thermodynamic entropy. On the gravity side this translates to the notion 
of isolated horizons as those for which entropy can be defined precisely. The framework of slowly evolving geometry 
is especially useful here, since it quantifies how far from equilibrium the black brane is by examining the validity of 
the first law of thermodynamics. This tool can be used as a probe of when and where the strongly coupled boundary 
quantum field theory is close to equilibrium. 



Note that in hydrodynamics due to the requirement of positive divergence of the entropy current, Ai and Bi contributions to the entropy 
current are not independent, but rather linked by the equation l|102|l . However such considerations were done under the assumption that 
the shear tensor vanishes locally, which is never the case in the Bjorken picture unless all the dissipative contributions are negligible. 
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Apart from some more or less obvious generalizations it would be interesting to explore these ideas in the context of 
equilibration of the boundary quantum field theory perturbed out of equilibrium by localized sources (in the spirit of 
12, 36, This application of the slowly evolving horizons framework is especially interesting because in the case 
of planar horizons considered here there can be widely separated regions, of which some are in local equilibrium while 
others are not. 

On the gravity side these ideas also raise some interesting questions. Throughout this paper we have often referred 
to various surfaces being "close" to a future outer trapping horizon (FOTH). In practice this has always meant close 
along surfaces of constant Eddington-Finkelstein coordinate v in terms of the radial coordinate r of the adapted 
coordinate system. From a purely geometric point of view this is not really satisfactory; notions of distance shouldn't 
be formulated in terms of an arbitrary coordinate system. Thus, more thought needs to be given to properly under- 
standing whether there is an invariant way to define this closeness. One such idea has been addressed in [28] where the 
authors considered the minimum proper time interval between a spacelike FOTH and a null event horizon. However 
the results are not entirely satisfactory and work remains to be done. Ideally one would like to obtain a general proof 
showing that in the vicinity of a slowly evolving horizon there would always be candidate event and alternate horizons 
which are arbitrarily close in a well-defined manner. How to do this is currently an open problem. 

This work also sheds some light on the long-standing discussion as to whether it is more "correct" to consider event or 
apparent horizons in physical situations. Instead of there being just two choices, we now propose that are many more 
and the uncertainty as to which one should considered actually reflects a physical ambiguity in the proper definition 
of entropy and other quantities such as energy as one moves away from equilibrium. Again these issues deserve further 
investigation. 
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APPENDIX A: DERIVATION OF THE CONSTRAINT EQUATIONS 

In this appendix we consider the derivations of the geometric constraints mentioned at the end of section IIIBI A 
more complete discussion of almost all of the results mentioned in this section can be found in j2o| . 

To begin we write V in the more general form 

which will help to display symmetries in the equations. In the main text A is set back to unity. 

Now, in finding the equations constraining the derivatives of the extrinsic geometry the key piece of information that 
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must be used is that la and Ua are always normal to the S\ . Thus 



qlCvh = =^ V^^da - -{da- ^a)C + Kyia and 

q^a^V^b^Q =^ V''Vbna= {da+LUa)A- Kvna, 



(A2) 
(A3) 



The importance of these two relations is that they can be used to replace covariant derivatives in the V° direction 
with derivatives in the tangent direction plus a part proportional to Ky: the V° component of the connection on the 
normal bundle. 

It is then possible to apply these relations to find expressions for the rate of change of components of the extrinsic 
geometry. For fl^^-j one proceeds in the following manner. First, 



Expanding the first term with ([3]) and then applying (jA2D and (IA3p gives 

(V^Vcr')Va4 = {d''A + w''A){eVaib) - id''C - Cj''C)[n''Vatb) 



(A4) 



(A5) 



while for the second term, the Riemann tensor can be used to commute the indices of the double derivative to obtain 

r'V^V.Va^ = q'^^'V^in^abdi" + ^a^ch) ■ (A6) 

Next, applying the definition of the Einstein tensor Gab — T^ab — ^T^9ab along with the Gauss relation 

ilqUltan^f^H = Rabc, + {ki'^ki:^ + k^l^) - [k^k^ + kt'Jk^l) , (A7) 

where Rabcd is Riemann tensor associated with the {n — l)-dimensional qab, the Riemann term of (jA6p becomes 

^'V^n^abJ'' = -AGabi^'i'' - CGabi'^n^ + Ir + 0^e)Oin) - k^ab^tl) ■ (A8) 
Meanwhile, again with the help of (|A2[) the other term is 

.(V).a6 



r'V^V^V.^ = -d'C + da{CCo'^) + nve(,)^k)^^,'kl^^ 

Combining these results there is some cancellation and so one gets 

CvOi^i) ^ Kv6(^- d^C + 2Cj''daC -c\\\Cjf -daCb" ^ R/2 + Gabl°"n!' -0(i)0(n) 

— A k^ab^ti) + Gabi'^C 



(A9) 



(AlO) 



where \\uj\\'^ = uiaCo"' ■ In this form the relation is independent of the spacetime dimension, however, if one decomposes 
the extrinsic curvature term on the second line with ([5]) the dimension appears 



CvO{i) = Kve(^i)-d^C + 2Cj''daC -C 



-A 



3l 



(All) 
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where and Hcr'-^-'lP = cr^^hCr'^iy Note in particular that the only quantity not directly defined on Sy is the gauge- 
dependent Ky. 



The rate of change of the inward expansion may be obtained simply by exchanging i and n and swapping A ^ —C 
in mi). Then 



anb 



n — 1 



(A12) 



A similar calculation using (|A2p and (jASp finds the derivatives of the expansions pulled back onto the Sx . First again 
using the Riemann tensor to commute derivatives one can show the Codazzi relations 



^:q[q'j''n,fgh = {da-Cja)ki'J -{db-Cub)ki'^ 



(A13) 



"-gt^^'^e/,. = {da + Cja)ki:^ ~ {db + C^b)k^:'} 



Now, the extrinsic curvatures can be decomposed into expansions and shears and while the Riemann tensor can be 
decomposed into Weyl and Ricci components using the well-known relation [18| : 



2 2 

Tlabcd = Cabcd + {ga[c'Ti-d\b + fff)[d'^c]a) ^ ^^^^ _ {ga[cgd]b) 



(A14) 



where Cabcd is the Weyl tensor, square brackets indicate anti-symmetrization and the reader should keep in mind that 
we are considering an {n + l)-dimensional spacetime. Then constracting with (f"^ we find 

1 



daO^n) = -e^n^C^a^idb+Qbhi"^^' 



q'.GbJ' - 2q'^Cbcdei'£''n^ and 



cnd^e 



n — 1 



1 



qlGbcn' + 2qlCbcden^tri 



n-l 



(A15) 
(A16) 



The final equation of interest for this paper is the rate of change of the angular momentum one- form up A. This 
comes from a direct expansion of q^C\>{ncVbi'') along with the familiar applications of the Riemann tensor to commute 
derivatives, (jA2p and (|A3|) to enforce the fact that the S\ fit together smoothly into a surface, and then breaking the 
Riemann tensor into Weyl and Ricci components. We find 



CyC^a = daKv - fc^ [d" A + [u" A] - fc^,'^ [d" C - u'C] 

■^^—Gbcir + Cn') - CbcdeV't'^n' 
_(n- 1) 



(A17) 



As shown in the main text, these equations lie at the root of black hole mechanics and dynamics. 
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APPENDIX B: BJORKEN FLOW AT SECOND ORDER 



The second order solution 



26 



27| for the coefficients in ((57)) reads 



A2 {v) 



b'2 iv) 



2tt^ log(w)A2 



1 A2 



27r2 log(7rA)A2 



3w4 3?;4 3^4 

TT^ log(2)A2 {v^ - Stt^A^) -2^;4 + ^^A^y^ + 4^3y^3^ ^ 2tt^A^ 

2 (3v4 - 27r3A3w - 37r4A4) Si 2 (w" + tt^A^) ,52 



+ 



9u6 

^4A4)tan-i(^) 



3w5 



37ru^A 
!log(i;)w^ 



! log(7rA)w2 {7v^ - Att'^A^v - 3^^^^^ tan"! (^) 



97r5A5 - Ovrw^A gyr^AS - Ottv^A OttA (w^ - 7r4i;2A4) 

4^2 A2 log(2) 2 {v^ + ttAv + n'A^) 2 (v^ - ttAv + n^A^) log + 
97r4uA4 - 9^5 9A(t; + ttA) (v'^ + tt^A^) w 97rA(7rA - v) {v^ + tt^A^) v 

4 (t;2 + ttAv + ^^A^) log + l) 2^ 25| _ 

97rA(u + ttA) (^2 + 7r2A2) t; ^ 3^2 9t;3 
11^6 + 227rAz;5 + 347r2A2i>4 + QOn^A^v^ + 437r4A4i;2 + 267r5A5i; + 1677*^ A^ 

9(f + 7rA)2 (z;2 + 7r2A2)%3 ^ 
4 (3u^ + 67rAw5 + 9tt^A'^v^ + Itt^A^v^ + 5t:'^A^v^ + 3t:^A^v + Tr^^A^) Si 



9(i; + 7rA)2 (1,2 + 7r2A2)^ i;3 



C2 (w) = d2 (v) - b2 (v) 
d2 (v) 



SI Si 


V — ttA 


6w2 3w2 


37rw2A 


log(^ + l) 




47r2A2 


6 



67r2A2 



67r2wA2 

21og(^A) S2 
37r2A2 w 



21og(t^) 
37r2A2 

1 

^ 127rA2 



(Bl) 



(B2) 



(B3) 
(B4) 



The expression for 62 given above can be integrated in terms of polylogarithmic functions, but the explicit form of 
the integral will not be needed in the sequel. The quantities Si (and S2.3,... in higher orders) are integration constants 
related to residual gauge symmetry r ^ r + R (f ) mentioned earlier. 
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